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We study the steady state of diffusion-limited coalescence, A + A ^ A, in the presence of a trap 
and with a background drift. In one dimension this model can be analyzed exactly through the 
method of inter-particle distribution functions (IPDF). Because of the irreversible trap the steady 
state of the system is a non-equilibrium state. An interesting phase transition, controlled by the 
drift away from the trap, takes place: from a non-trivial steady state, when the drift is weak, to 
a trivial steady state (the vacuum), as the drift increases beyond some critical point. Surprisingly, 
regardless of the drift strength, the computed IPDF resembles that of an homogeneous equilibrium 
system, without the trap. We suggest that this is due to "shielding": the particle nearest to the trap 
shields the remaining particles from the effects of the trap. Finally, we compare the exact solution 
to that of a reaction-diffusion equation, and we determine the optimal values of the appropriate rate 
coefficients. 



I. INTRODUCTION 

Non-equilibrium kinetics of diffusion-limited reactions 
has been the subject of much recent interest [0-0. 
In contrast to equilibrium systems — which are best 
analyzed with standard thermodynamics — or reaction- 
limited processes — whose kinetics is well described by 
classical rate equations — there is no general ap- 
proach to non-equilibrium, diffusion-limited reactions. 

In this paper we study a diffusion-limited coalescence 
process in one dimension: A + A ^ A, which yields it- 
self to exact analysis |lC|-|2l|] . We show that a solution is 
possible even in the presence of background drift, and in 
the presence of a trap; when in the long-time asymptotic 
limit the stationary state of the system is a nonequilib- 
rium state. The drift field controls an interesting phase 
transition in the steady state of the system: The size of 
the depletion zone next to the trap increases as a func- 
tion of the drift away from the trap. Beyond a certain 
critical drift strength, the only possible stationary state 
is the vacuum (an empty lattice, with no particles). 

Our analysis of the steady state yields the surprising 
result that the distance between a given particle and the 
nearest particle to its right (the "forward" IPDF) is ex- 
actly the same as in equilibrium, even in the presence of 
a trap, and independently of the drift! We suggest that 
this is due to a "shielding" effect: the particle nearest to 
the trap shields the rest of the particles from the effect 
of the trap. The trap and the drift influence only the 
distance between the trap and the nearest particle to the 
trap — all other particles remain unaffected. 

Finally, we use the coalescence model with a trap to 
study the applicability of reaction-diffusion equations. 
Reaction-diffusion equations are the most common ap- 



proximation method in the study of diffusion-limited pro- 
cesses. Our system is a rare example where (in the 
absence of drift) the stationary solution of the relevant 
reaction-diffusion equation may be found in closed form. 
We explore in what ways the reaction-diffusion equation 
is approximate, by comparing its results to the exact so- 
lution of the IPDF method. The comparison also al- 
lows us to determine the optimal rate constants of the 
reaction-diffusion equation in a straightforward manner. 
Normally, this feat requires a renormalization group anal- 
ysis. 

The rest of this paper is organized as follows. In Sec- 
tion [n| we present a lattice model of diffusion-limited re- 
versible coalescence, along with the exact method of anal- 
ysis; the method of Empty Intervals, also known as the 
method of Inter-Particle Distribution Functions (IPDF). 
The stationary equilibrium state in a homogeneous in- 
finite system, needed for comparison with the trap, is 



analyzed and summarized in Section III. In Section IV, 



a trap is introduced and the resulting stationary state 
is explored. The exact solution includes a description of 
the dynamical phase transition in the stationary kinetics 
of the system, brought about by the background drift. 
The explanation to the surprising result that the IPDF 
of the system with a trap is homogeneous (and exactly 
the same as in equilibrium) is presented in Section [v|. 
Last, in Section |Vl| we compare the exact solution for 
the case of non-biased diffusion (zero drift) to that of the 
corresponding reaction-diffusion equation. Our goal is to 
show in what ways the latter is an approximation, and to 
devise strategies to find out appropriate rate coefficients. 
We conc lude with a discussion and open questions in Sec- 
tion IvlI. 
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II. REVERSIBLE COALESCENCE 

Our model |l^| is defined on a one-dimensional lat- 
tice of lattice spacing a. Each site is in one of two states: 
occupied by a particle A (•), or empty (o). Particles 
hop randomly to the nearest neighbor site on their right 
(left), at rate D/a 2 + u/2a {D/a 2 - u/2a). Thus, in the 
continuum limit of a — ► 0, the particles undergo diffusion 
with a diffusion constant D, and with a uniform back- 
ground drift velocity u. A particle may give birth to an 
additional particle, into a nearest neighbor site, at rate 
v/a (on either side of the particle) p2|. If hopping or 
birth occurs into a site which is already occupied, the 
target site remains occupied. The last rule means that 
coalescence, A + A — > A, takes place immediately upon 
encounter of any two particles. Thus, the system models 
the diffusion-limited reaction process 



A + A^A 



(1) 



The dynamical rules of the model are illustrated in Fig. 1. 
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Figure 1. Reaction rules: (a) diffusion; (b) birth; and co- 
alescence, (c) following diffusion, and (d) following a birth 
event. The dotted lines in (a) and (b) indicate alternative 
target sites. In (a), the rate of the two possibilities differ due 
to the drift. 

An exact treatment of the problem is possible through 
the method of Empty Intervals, known also as the 
method of Inter-Particle Distribution Functions (IPDF) . 
The key concept is E n>m (t) — the probability that sites 
n, n + 1, • ■ ■ , m are empty at time t. The probability that 
site n is occupied is 

Prob(site n is occupied) = Prob( • ) = 1 — E n , n . (2) 

The event that sites n through to are empty (prob. E nim ) 
consists of two cases: site to + 1 is also empty (prob. 
En.m+i), or it is occupied. Thus the probability that 
sites n through m are empty, but site to + 1 is occupied 
is 



Prob(o • • ■ o •) = E n<m - E n>m+ i , 
and, likewise, 

Prob(» o • • • ™) = E n;m - B„-i,m • (4) 



(3) 



With this in mind, one can write down a rate equation 
for the evolution of the empty interval probabilities: 
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(•5) 



For example, the first term on the r.h.s. of Eq. (||) ac- 
counts for the increase in E n ^ m when the particle at the 
right edge of o ■ • • o • hops to the right and the sites 
n, . . . , to become empty; the second term denotes the de- 
crease in E n>m when a particle at to + 1 hops to the left 
into the empty interval n, . . . , m, and so on. 

Eq. (^) is valid for to > n. The special case of to = n 
corresponds to E n , n — the probability that site n is empty. 
It is described by the equation 



dt x a 2 ' 2a 
D u 
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(6) 



Comparison with Eq. (||) yields the boundary condition: 

E n ,n-1 = 1 ■ (7) 

The fact that the {-E„ jm } represent probabilities implies 
the additional condition that 



En,rn ^ ■ 

Finally, if the system is not empty then 



lim E„ 







(8) 



(9) 



In many applications, it is simpler to pass to the con- 
tinuum limit. We write x — na and y — ma, and replace 
E n ,m{t) with E(x,y;t). Letting a — * 0, Eq. (||) becomes 

dE . d 2 d 2 . .8E , ,dE 

(10) 

with the boundary conditions, 

E(x,x;t) = l, (11) 
E(x,y;t)>0, (12) 
lim E(x,y;t) = . (13) 
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The concentration of particles is obtained from the 
continuum limit of Eq. (g) , together with Eq. ([ll]) : 



dE{x,y\t) ^ 

C(X,l) — \y = x ■ 



(14) 



The conditional joint probability for having particles at 
n and m but none in between, is (Fig. 2) 



P n ,m(t) = Prob(2 o • • ■ o •) 

= ^n+l,m — 1 -E'n+l.m ^n,m— 1 "I" E nm , (15) 



which in the continuum limit becomes 

d 2 £(x,2/;i) 



P(z,y;t) 



9x dy 



(16) 



From P we can compute the IPDF. The "forward" IPDF, 
i.e., the probability that given a particle at x the next 
nearest particle to its right is at y, is 



P{x,y;t) = c(x,t) 1 P(x,y;t) 



(17) 



Likewise, the "backward" IPDF — the probability that 
given a particle at y the next nearest particle to its left 
is at x — is 



q(x,y;t) = c(y,t) 1 P(x 1 y;t) 



(18) 
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Figure 2. Discrete IPDF: E n +i,m-i consists of four events, 
listed above the "total" line. In this way one can find P n ,m. — 
the only remaining unknown. 



III. THE EQUILIBRIUM STATE 



We now derive the steady state solution of Eq. (10), 
with the boundary conditions (O) - ( 13). It proves useful 



to change variables: 77 = x + y, and £ = y — x. £ repre- 
sents the length of the intervals, and r\ represents their 
location (the center of the interval is at 77 / 2) . In these 
new variables, the stationary limit of Eq. (|lfj) becomes 



d 2 d 2 



dr] 



d_ 



(19) 



where we have used the shorthand notation u = u/D 
and 7 = v/D. 



Because the background drift is homogeneous, we ex- 
pect the solution to be translationally invariant, that is, 
the probability of empty intervals should depend only on 
their length but not upon their location: E = E(£). E 
then satisfies a simple ODE; d 2 E/d£ 2 — jdE/d^ = 0, 
whose solution (in view of the boundary conditions) is 



E eq = < 



This corresponds to the particle concentration (Eq. |lj) 

(21) 
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The subscript "eq" emphasizes the fact that here the sta- 
tionary solution is an equilibrium solution. This is to be 
expected because all the microscopic rules involved are 
reversible. Notice that the equilibrium solution is inde- 
pendent of the drift. The effect of the drift is to im- 
part the whole system a uniform background velocity, 
u, which in fact may be eliminated through a Galilean 
transformation. 

We can test that the distribution of particles is in- 
deed an equilibrium distribution, by working backwards. 
At equilibrium the particles ought to be distributed ran- 
domly, independently from each other — a state of maxi- 
mum entropy. Suppose that the particle concentration is 
7, then the probability that there are no particles in an 
infinitesimal interval of length A£ is (1 — 7A£). Since the 
particles are uncorrelated, the probability that a finite in- 
terval of length £ = y — x is empty is then (1 — 7A£)£/ A ^ . 
In the limit A£ — > 0, one recovers the equilibrium empty 
interval probability, Eq. (|20|). 

Finally, we can work out the IPDF. Using Eqs. @ - 
(|l8|), (pb[), and (pl|), we find the Poisson distribution: 



Peg(x,y) = q eq (x,y) = -ye 



— ryp--yiV- x ) 



(22) 



Again, this can be obtained, independently, from an anal- 
ysis of uncorrelated, randomly distributed particles. 



IV. STEADY STATE WITH A TRAP 



The equilibrium solution, Eq. (|20|), is not the only solu- 
tion to Eq. ( p9| ) . Another trivial solution is the "vacuum" 
state; E(x, y) = 1. It represents a completely empty sys- 
tem. But the solution ( pC| ) is stable, while the vacuum 
state is not. In fact, when the initial state of the system 
is a mixture of the two phases: c(x,t = Q) = for x < 0, 
and c(x, t — 0)— c eq for x > 0, say, then the stable phase 
invades the unstable phase. In the absence of drift, the 
front between the two phases propagates at the constant 
speed v (the birth rate) |2^] . This system has been stud- 
ied as a "noisy" analogue of the mean-field problem of 
Fisher waves |2J]. Here we wish to study another inho- 
mogeneous situation, where there is a perfectly absorbing 
trap at the origin, instead of the initial empty half-space. 

As opposed to the Fisher problem, in the case of a 
trap the system arrives at a stationary state. This state 
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must be a non-equilibrium state, because the trap is ir- 
reversible (particles trapped in can't ever come out). As 
the trap depletes its immediate neighborhood, the deple- 
tion zone is continually replenished by a stream of par- 
ticles from the stable phase which rush in at speed v, as 
in the Fisher wave problem. Suppose now that the par- 
ticles are subject to a drift u away from the trap, then 
an interesting transition would take place as u is made 
larger than v: the depletion zone would grow faster than 
it could be replenished, and the stationary state would 
never be achieved! We wish to study this transition, 
and how the depletion zone is affected by the drift when 
u < v. 

To derive the appropriate boundary condition, we turn 
back to the discrete representation. In the presence of a 
perfect trap at n = 0, Eq. (||) is limited to < n < m. 
The special equation for n = is 

dEo.„ 
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(23) 



Comparison to Eq. 
dition 



yields the discrete boundary con- 



which in the continuum limit becomes 
dE(x,y;t). 



(24) 
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. 



In addition, the boundary condition ( |l3| ) is now replaced 
by 



lim E(0,y;t) = 

y— too 



(26) 



We are looking then for the solution to Eq. (19|), con- 
fined to the wedge (0 < x < y), and which satisfies 
the boundary conditions Q), @, @, and @. Fol- 



lowing the method of separation of variables, assume 

E{r},0 = Ea^/oWfcK). tncn fa and g a satisfy the 
ODEs: 



d 2 f a 
dr] 2 
d 2 g a 

de 
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dfa 
! 
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dg a 
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(27) 
(28) 



We find the solution by inspection: A good guess is that 
a 2 = would be part of it, since far away from the trap 
one expects convergence to the equilibrium state. How- 
ever, /ogo fails to satisfy the boundary condition due to 
the trap, Eq. (p^). Looking for the simplest possible way 



to mend this and to satisfy the remaining boundary con- 
ditions, we find that it is sufficient to consider just one 



additional component: 4a 2 
turns out to be 
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The final solution 



E s (x,y) = e-Tto"*) + l e -TV( e "W - e^) 



(29) 



Far away from the trap, as x,y-^ oo, this converges to 
the equilibrium result of Eq. ( |20| ) . 

From (|l4|), we obtain the stationary concentration pro- 
file: 

c s (x) = 1 [l-e-^-^ x ] . (30) 

There is a depletion zone of size (7 — w) _1 = D/(v — u) 
near the trap, and the concentration grows asymptoti- 
cally to c eq — 7 as x — y 00. As the drift velocity u grows 
and approaches v, the depletion zone becomes larger and 
larger. In the limit u —> v the depletion zone is infinite: 
from (^9| ) we see that then E s — » 1, i.e., the stationary 
state is the vacuum (but it takes an infinite time to get 
there). 

The IPDF is surprising. From (jl^) we get the condi- 
tional joint probability 

P s (x, y) = ! 2 e-^y- x \l - e -(7—)«] 

and so the forward IPDF is 

-7? 



(31) 



Ps(x,y) = je~ 



£, = y-x 



(32) 



regardless of the drift. The notation chosen here empha- 
sizes the fact that p s is translationally invariant, in spite 
of the trap at the origin. In fact, the forward IPDF is 
exactly the same as in the equilibrium state (Eq. p2| )! 
The backward IPDF is 



(25) q s (x,y)=-f- 



-vSp-h-^y 



1 



-(y-u)y 



(33) 



Notice this time the dependence on the position of the 
given particle, y, as well as the dependence on the drift, 
and that the backward and forward IPDFs are not equal. 
Moreover, the backward IPDF does not normalize prop- 
erly. The reason for that is that there is a finite chance 
that there are no particles between the particle at y and 
the trap; i.e., the particle at y might be the nearest par- 
ticle to the trap. The probability density for the distance 
of the nearest particle to the trap can be computed from 
the continuum limit of Eq. (j}|); po(y) — ~dE/dy\ x= a, or 



Po(y) 



7 [e ^ 



■j)y 



-72/1 



(34) 



With this understanding, the proper normalization con- 
dition is 



q(x,y)dx + c(y) 1 p (y) = 1 



(35) 



which is indeed met! 

As a last remark, notice that all of the above results are 
also valid when the drift is toward the trap (u < 0). In 
this case the depletion zone shrinks as the drift's strength 
increases, and it vanishes in the limit of infinite drift. 
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V. SHIELDING 

At first sight there seems to be a contradiction be- 
tween the forward IPDF (Eq. |3^) and the concentration 
profile (Eq. ^). How could it be that the distribution 
of distances between nearest particles is independent of 
location, while particles are sparser near the trap!? - 
The answer is that the forward IPDF is a conditional 
probability — dependent on the presence of the first parti- 
cle at x. Thus, the unexpectedly small distance between 
particles in a sparse region is compensated by the rare 
likelihood of finding such pairs of nearest particles in the 
first place! Nevertheless, one can't help wondering what 
kind of distribution would explain the results of the pre- 
vious section. The answer turns out to be surprisingly 
simple. 

Suppose that the particles to the right of the trap are 
distributed just as in equilibrium, with the exception of 
the nearest particle to the trap. Assume further that the 
nearest particle to the trap is at distance z, with proba- 
bility Po{z) (the same probability density as in Eq. 34). 
In this view, the effect of the trap is merely to create 
an empty gap between itself and the first particle — a 
gap which increases as the drift u approaches the crit- 
ical value v. [It follows from Eq. (f34|) that the average 
distance to the first particle is (2v — u)/{v — u).] Another 
way of looking at the suggested distribution is as if the 
nearest particle to the trap shields all the other particles 
from the effects of the trap! The proposed distribution 
is illustrated in Fig. 3. 

We now show that this model reproduces the results 
All we need to do is to derive the empty 



of Section IV 



interval probability E[x, y) implied by the model, since 
everything else follows from it. E(x,y) can attain one of 
three values — depending on the location of the endpoints 
of the (x, y) interval relative to the location of the nearest 
particle to the trap, z: 



1 x < y < z 

E(x, y) = i x < z < y 



-~i(y-x) 



(36) 



z < x < y 



Therefore, 



E{x,y)= p (z)dz + e-^y-^ Po (z)dz. (37) 

Jy JO 

The result is exactly as E 8 of Eq. (p9|)! 



Trap 



nearest 
particle 



Figure 3. Schematic illustration of the shielding effect: The 
particles in the shaded area are distributed randomly and 
independently from each other, as in equilibrium. The gap 
z between the particles and the trap follows the probability 
density distribution po(z). 



So far, we have merely demonstrated that "shielding" 
is consistent with the exact results. A complete proof re- 
quires showing that the full hierarchy of n-point density- 
density correlation functions are the same as postulated 
by the shielding model. The problem of n-point cor- 
relation functions can be formulated exactly within the 
framework of the IPDF method . With this method, 
it is rather straightforward to prove the shielding ef- 
fect p5l. 



VI. REACTION-DIFFUSION EQUATION 

Reaction-diffusion equations are one of the most valu- 
able and widely used approaches to kinetics of reaction 
processes. The method posits the existence of a meso- 
scopic length scale within which the system is homoge- 
neous and well mixed, and where the reaction rates can 
be accounted for as in classical rate equations. At longer 
length scales, variations in the concentration, c(x, t), give 
rise to diffusion, which is modeled by a simple diffusion 
term. For example, the reaction-diffusion equation that 
applies to our model, in the absence of drift, is 



at ox z 
together with the boundary condition: 

c(0,t) = 0, 



(38) 



(39) 



which is imposed by the trap at the origin. D' is the 
effective diffusion constant, and fci and k 2 represent the 
effective rates of birth and coalescence of particles, re- 
spectively. 

In spite of the popularity of reaction-diffusion equa- 
tions it is not often appreciated that they are mere ide- 
alizations, yielding only approximate results. The reason 
for this misconception is twofold: (1) few exact solutions 
of reaction-diffusion models are available for comparison, 
and (2) the reaction-diffusion equations themselves are 
usually complicated enough that they cannot be solved 
in closed form, and so their consequences are not always 
clear. Diffusion-limited coalescence with a trap is a rare 
example where the stationary state is known exactly, and 
where also the corresponding reaction-diffusion equation 
may be solved exactly, at least in the case of zero drift 
p6[ . We wish to compare the two solutions. 

The exact steady state for zero drift, needed for the 
comparison, may be obtained from the results of Sec- 
|[V|. Taking the limit w -*• in Eqs. (p|) and 



tion 



we find 



and 



E,(x, y) = e - l(y - x) + j(y - x) e - ry 



c s (x) = 7(1 



(40) 



(41) 
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Regarding the reaction-diffusion equation, Eq. (38), 
our first concern is to identify the constants D', k%, and 
ki- We note first that 



VII. DISCUSSION 



D' = D , (42) 

since in the absence of reactions (when k\,k2 = 0) the 
particles perform simple diffusion, characterized by the 
diffusion constant D — the same D as in the hopping rate 
D/a 2 of the microscopic rules. Furthermore, in the con- 
tinuum limit the model's only physical parameters rele- 
vant to its kinetics are the diffusion constant D (of di- 
mension L 2 /T; L = length, and T = time) and the birth 
rate v (dimension L/T). Therefore, the only way to pro- 
duce the required dimensions of k% (1/T) and ki (L/T) 
is by having k\ ~ v 2 / D and ki ~ v. Finally, consider the 
stationary solution of Eq. ( |38| ) in an infinite system with- 
out the trap: c eq = k\/ki- To conform with the exact 
solution of jj~t| ) we must have kx/kq, = v/D. Thus, 

ki = f— ; k 2 =0v, (43) 

where (3 is a constant. 

The stationary solution to Eq. (|3^) with the trap — the 
boundary condition (|||) — is 

= I t an h 2 U[*L X + tanh- 1 - - , (44) 

Coo 2 IV 3D V 3 / 2 ' v ' 

where Coo = k\/ki = v/D is the concentration of parti- 
cles infinitely far away from the trap. The concentration 
profile described by Eq. (H) looks similar to the exact 
result of Eq. (f4l|). One could now use different criteria 
to determine the value of the fitting parameter (3. De- 
manding the correct asymptotic behavior far away from 
the trap; lim^^oo ln[l — c s (x)/c eq ]/x = —v/D, we get 
(3 = 3/4. On the other hand, if we require the correct 
behavior close to the trap; (dc s /dx) x= o = (v/D) 2 , we get 
(3 = 9/4. Clearly, it is impossible to find a value of (3 that 
reproduces the short range behavior and the long range 
behavior simultaneously. In Fig. 4 we show the results of 
a least square fit in the range < (v/D)x < 5, which is 
achieved with (3 — 1.27 . 
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D 

Figure 4. Comparison of the exact analytic concentration 
profile (solid line) to the prediction of the reaction-diffusion 
equation (broken line). Shown is the result of a best fit, 
achieved with f3 = 1.27 . 



We have solved the one-dimensional problem of 
diffusion-limited reversible coalescence with a trap, and 
with a background drift, exactly. The system exhibits a 
dynamical phase transition in its long-time asymptotic 
steady state, which is controlled by the strength of the 
drift away from the trap. As the drift increases the de- 
pletion zone near the trap grows larger, until when the 
drift crosses a critical threshold and the depletion zone 
becomes infinitely wide; the system is then empty. 

A very intriguing prediction of the exact solution is 
the distribution of distances between any given parti- 
cle and its nearest neighbor to the right (the forward 
IPDF), which turns out to be exactly the same as in 
equilibrium — as in the system without the trap — and in- 
dependent of the drift! This is explained by a "shielding" 
effect: the nearest particle to the trap effectively shields 
the rest of the particles from the trap's influence. The 
trap and background drift affect only the width of the 
gap between the trap and the nearest particle. 

We have also contrasted the exact solution, in the ab- 
sence of drift, with the alternative, traditional approach 
of reaction-diffusion equations. Our model is a non- 
trivial example where both the true kinetics, as well as 
the solution to the model's reaction-diffusion equation, 
are known exactly. Comparing the two solutions we were 
able to relate the effective rates of the reaction-diffusion 
equation to the microscopic rates of the underlying pro- 
cess, without appealing to renormalization group tech- 
niques. 

There remain several interesting open questions. The 
exact solution of Section |^ is basically an inspired guess. 
It would be useful to develop a more formal derivation 
method — one that would allow further analysis of differ- 
ent generalizations of the problem at hand. Other open 
problems include the question of transient behavior. For 
example, beyond the transition point, when the drift ex- 
ceeds the critical value of u c = v, it takes an infinite 
time for the system to achieve the empty steady state. 
What is the exact time dependence at different points 
of the phase diagram? — What is the mean-field critical 
dimension beyond which the reaction-diffusion equation 
becomes exact? Previous work on a closely related prob- 
lem (Fisher waves [|i|) has suggested that d c = 3, but it 
was based on computer simulations of an "infinite" sys- 
tem and the results were controversial. The bounding of 
our present model by the trap at the origin may offer a 
better alternative for numerical studies. — What is the 
prediction of the reaction-diffusion equation formulation 
when drift is included? How well is the dynamical phase 
transition captured in this approximation? These and 
similar issues will be the subject of future research. 
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